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Abstract 

The Neotropical butterfly Dryas iulia has been collected from several locations in Thailand and Malaysia since 2007, and has 
been observed breeding in the wild, using introduced Passiflora foetida as a larval host plant. The butterfly is bred by a 
butterfly house in Phuket, Thailand, for release at weddings and Buddhist ceremonies, and we hypothesized that this 
butterfly house was the source of wild, Thai individuals. We compared wing patterns and COI barcodes from two, wild Thai 
populations with individuals obtained from this butterfly house. All Thai individuals resemble the subspecies D. iulia 
modesta, and barcodes from wild and captive Thai specimens were identical. This unique. Thai barcode was not found in 
any of the 30 specimens sampled from the wild in the species' native range, but is most similar to specimens from Costa 
Rica, where many exporting butterfly farms are located. These data implicate the butterfly house as the source of Thailand's 
wild D. iulia populations, which are currently so widespread that eradication efforts are unlikely to be successful. 
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introduction 

The introduction of exotic species to novel habitats is one of the 
most significant threats to biodiversity conservation. Introduced 
plants can become invasive, replacing natural vegetation. Intro- 
duced predators can consume indigenous prey that lack suitable 
defenses, and introduced insect herbivores can become plant pests, 
causing damage to native plants as well as crops [1]. Extensive 
import laws and quarantine procedures exist in nearly every 
country to curtail unintentional introduction of pestiferous insects, 
which could potentially "hitchhike" on imported plants or 
agricultural produce. Despite the tremendous effort spent trying 
to prevent the spread of insects between countries, introductions of 
insect pests are common [1,2]. 

Once a newly introduced insect has been discovered, it may not 
possible to determine how the introduction occurred or whether 
the same species was introduced multiple times. For example, 
Eastwood and coUeagxies [3] used DNA barcodes to demonstrate 
that all sampled Dominican Papilio demoleus shared a single 



barcode also found throughout most of Southeast Asia. Since the 
particular barcode haplotype found in the Dominican Republic is 
widespread in the Oriental Region, it was not possible to 
determine the precise location of the source population using 
DNA barcodes or whether the species was introduced more than 
once [3]. Determining the provenance of this introduction is 
important because the Southeast Asian lineage of this species is 
frequently a pest of Citrus, whereas the lineage from Australia and 
New Guinea is not [4,5] . In addition to assessing the potential crop 
damage an introduced insect species may cause, knowledge of a 
species' home range might also be useful for identifying suitable 
parasitoid species for biological control. 

Beginning in 2007, several independent observers recorded 
specimens of the Julia butterfly, Dryas iulia (Fabricius, 1775), at 
several locations in Thailand and Malaysia (Fig. 1) including 
Samui Island (Surat Thani province; Les Day, pers. comm.), 
Phuket Island (Phuket province; Sin Khoon Khew, pers. comm.), 
Tioman Island, Malaysia [6], and Phi Phi Don Island, Thailand 
(Krabi province; D]h,pe7S. ohs.). Kuppers [7] reported the species 
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Figure 1. Collection localities of Dryas iulia butterflies on the 
Thai-Malay peninsula (unpublished data) [6,7]. Orange butterfly 
symbols indicate localities from which we sampled specimens for this 
study; black symbols indicate unsampled localities from which the 
species has been recorded. The image illustrates the wings of specimen 
DL-08-T033 caught on Phi Phi Don Island, Krabi Province, Thailand. 
doi:1 0.1 371 /journal.pone.01 04076.g001 

from the Thai provinces of Nakhon Si Tliammarat, Phang Nga, 
and Chumphon, and suggested that the species miglit have 
escaped from a butterfly house on Phuket. 

Dijas iulia is native to the Americas, where thirteen subspecies 
are found in the southern USA, Central America, the Caribbean, 
and northern South America [8,9]. To identify the Thai specimens 
to subspecies (Fig. 1), we compared wing patterns of this material 
to published photographs of all subspecies [9-1 1]. Specimens from 
Thailand resemble the subspecies Dryas iulia modesta, found in 
Texas, Mexico, Central America, and the Pacific coast of South 
America to Ecuador [8]. In their native range, larvae of Z). iulia 
feed on a variety of different Passiflora species (Passifloraceae). 
Plants in this genus are typically vines or lianas, with more than 
100 species in the New World tropics, and about 20 in tropical 
Southeast Asia, Australia, and New Zealand [12]. The second 
edition of Butterflies of Thailand [13] now lists the species as being 
part of the country's fauna, noting its presence in Chumpon, 
Nakhon Si Thammarat, and Surat Thani provinces, which are all 
in the southern peninsula. In Thailand, larvae have been found 
feeding on Passiflora foetida, and adults frequently nectar on 
Lantana camara (Les Day, pers. comm.); both of these plants are 
invasive species native to the Americas [14,15]. 

The Phuket Butterfly Garden (PBG; 71/6 Moo 5, Soi Paneung, 
Yaowarat Road, Rassada Rd., Phuket City) has been open since 
1990 in the center of Phuket Island, one of the country's most 
visited tourist areas. In addition to maintaining a butterfly 
vivarium, the company sells live butterflies for release at weddings 
(phuketbutterfly.com/weddrng.php, accessed March 2014) and 
provides butterflies for mass pubhc release [16,17]. Releasing 
butterflies at weddings is a relatively new custom practiced around 
the world. Instead of throwing rice or birdseed at newlyweds as 
they leave the wedding ceremony, celebrants release live butterflies 
from an envelope or cage so the couple departs in a swarm of live 
insects [18,19]. After noticing that the PBG website (phuketbutter- 
fly.com, accessed May 2008) showed pictures of D. iulia 
butterflies, a Thai colleague visited the facility at our request in 
June 2008. He found D. iulia flying in the vivarium, confirmed 
with staff that D. iulia could be purchased for release at weddings, 
confirmed that specimens could be shipped to the resort islands of 



Samui and Phi Phi Don (where D. iulia has already been observed 
in the wUd), and obtained nine fresh specimens of this species. We 
subsequendy froze the specimens for genetic work. 

In the present study, we address two questions: 1) Did wild Thai 
populations of Dryas iulia originate from livestock at the Phuket 
Butterfly Garden (PBG)? and 2) From where in its natural range 
did PBG animals originate? To answer these questions, we 
sequenced the barcodiiig section of the mitochondrial cytochrome c 
oxidase subunit I (COI) gene from wild-caught specimens in 
Thailand and specimens obtained from PBG. COI evolves rapidly, 
is easily amplified and sequenced with highly conserved primers, 
and is therefore a good marker for assessing maternal relatedness 
and, potentially, species membership [20]. If wild Thai D. iulia 
were naturalized after introduction from PBG, then DNA 
barcodes from wild-caught specimens would be similar or identical 
to barcodes from specimens obtained from PBG. Genetic 
differences between wild and PBG-derived specimens would 
suggest that PBG is not the source of Thailand's naturalized D. 
iulia. However, shared barcode sequences might also result if 
different populations of the species do not vary at this locus. 
Therefore, we compared these Thai sequences to barcodes from 
D. iulia modesta specimens sampled throughout the species' native 
range, including sequences from GenBank and from D. iulia 
hispaniola specimens wild-caught in the Dominican Republic. We 
suspect that PBG stock originated in Costa Rica, as many 
Neotropical butterfly farms are found here (Michael Boppre, pers. 
comm.) [18]. 

Materials and Methods 

Specimen acquisition 

Butterfly specimens were caught with an aerial net in the field. 
Each specimen's wings were removed from its body. Wings were 
stored in glassine envelopes and bodies were placed in vials of 
100% ethanol and frozen. AH specimens collected for this study 
are vouchered in the DNA and Tissues Collection of the Museum 
of Comparative Zoology, Harvard University (Table 1). Permis- 
sion to conduct research in Thailand was granted by the National 
Research Council of Thailand. Permission to export specimens 
was granted by the CITES Office of the Department of National 
Parks, Wildlife and Plant Conservation. Permission to conduct 
research in the Dominican Republic was granted by the Ministerio 
de Medio Ambiente y Recursos Naturales. Permission to export 
specimens was granted by the Secretaria de Estado de Agricultura, 
Departamento de Vida Silvestre, Santo Domingo. All permits in 
the Dominican Republic were arranged by Kelvin A. Guerrero 
(kguerrero.net). 

DNA sequencing 

Specimens were obtained from colleagues, from the Phuket 
Butterfly Garden (see Introduction), and from field collection in 
Thailand and the Dominican Republic. DNA was extracted from 
single butterfly legs using a QIAGEN DNEasy Blood & Tissue Kit. 
After addition of the tissue lysis buffer, insect legs were ground 
mechanically in microcentrifuge tubes using disposable pestles. 
This step was added to further break down the chitin exoskeleton 
and thereby maximize the surface area of tissues exposed to the 
lysis mixture. Subsequenfly, proteinase-K was added and the 
manufacturer's protocol was resumed. 

A 658 bp fragment of the mitochondrial gene cytochrome c 
oxidase subunit I (COI) was amplified from whole genomic 
extracts using the diverse metazoan invertebrate primer pair 
LCO1490 (5'- TGTAAAACGACGGCCAGTGGTCAA- 
CAAATCATAAAGATATTGG-3') and HC02198 (5'-CAG- 
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Table 1. Specimen Information for sequences included In this analysis. 



Species 


Location: 
Voucher # 


Coll. Date 


Coordinates 


Collection Locality 


GenBank 


Dryas iulia hispaniola 


MCZ: NC-n-J036 


22-Apr-l 1 


19.0676 N 70.8639 W 


Dominican Republic: La Vega, La Cienaga 


KJ496350 


Dryas iulia hispaniola 


MCZ: NC-11-J053 


23-Apr-l 1 


19.0969 N 70.6186 W 


Dominican Republic: La Vega, Salto Baiguate 


KJ496351 


Dryas iulia modesta 


DHJ: 04-SRNP-56202 


l-Jun-08 




Costa Rica 


GU157122 


Dryas iulia modesta 


DHJ: 04-SRNP-56042 






Costa Rica 


GU157127 


Dryas iulia modesta 


DHJ: 05-SRNP-58644 






Costa Rica 


GU157128 


Dryas iulia modesta 


DHJ: 05-SRNP-3241 1 






Costa Rica 


GU157129 


Dryas iulia modesta 


DHJ: 05-SRNP-45212 






Costa Rica 


GU157130 


Dryas iulia modesta 


DHJ: 05-SRNP-32379 






Costa Rica 


GU157131 


Dryas iulia modesta 


DHJ: 08-SRNP-55873 


5-May-08 


10.763 N 85.413 W 


Costa Rica 


GU666775 


Dryas iulia modesta 


DHJ: 08-SRNP-72539 


lO-Oct-08 


10.997 N 85.397 W 


Costa Rica 


GU666782 


Dryas iulia modesta 


L)HJ: U/-bKNr-4D9DD 






Costa Rica 


'J CO no 


Dryas iulia modesta 


DHJ: 03-SRNP-11685 


27-Jul-03 


10.9066 N 85.2878 W 


Costa Rica: Alajuela, Area de Conservaclon Guanacaste, 
Rincon Rainforest, Sendero Juntas 


GU333908 


Dryas iulia modesta 


DHJ: 03-SRNP-16486 


18-Jul-03 


11.0182 N 85.4502 W 


Costa Rica: Guanacaste, Area de Conservaclon Guanacaste, 
Del Oro, Tangelo 


GU333906 


Dryas iulia modesta 


DHJ: 04-SRNP-49035 


7-Oct-04 


10.89 N 85.48 W 


Costa Rica: Guanacaste, Area de Conservaclon Guanacaste, 
Sector Cacao, Quebrada Otilio 


GU157124 


Dryas iulia modesta 


DHJ: U4-bKNr-4/oUD 


22-Aug-u4 


1 n on M oc AO lA/ 
1 u.o9 N OD.4o W 


Costa Rica: Guanacaste, Area de Conservaclon Guanacaste, 
Sector Cacao, Quebrada Otilio 


GUId/126 


Dryas iulia modesta 


DHJ: 04-SRNP-49966 


l-Dec-04 


10.886 N 85.482 W 


Costa Rica: Guanacaste, Area de Conservaclon Guanacaste, 
Sector Cacao, Sendero Guayabal 


GU157123 


Dryas iulia modesta 


MAL-02661 


26-Aug-99 


18.731 N 89.395 W 


Mexico: Campeche, Calakmul, Reserva de la Biosfera de 
Calakmul, 24 Km N X-PujII, Entrada a 'El Papagayo' 


GU659625 


Dryas iulia modesta 


MAL-02660 


14-Aug-02 


17.972 N 89.358 W 


Mexico: Campeche, Calakmul, Reserva de la Biosfera de 
Calakmul, Zona K, Dos Naciones 


GU659624 


Dryas iulia modesta 


MLL-00747 


21-NOV-06 


1 8.68 N 89.40 W 


Mexico: Campeche, Calakmul, Zoh Laguna 


JN201279 


Dryas iulia modesta 


MLL-01591 


29-Mar-07 


18.68 N 89.39 W 


Mexico: Campeche, Calakmul, Zoh Laguna 


JN201280 


Dryas iulia modesta 


MAL-02658 


14-May-04 


1 8.484 N 89.899 W 


Mexico: Campeche, Ejido Conhuas, Campamento Yax'che, 
Guardaraya Sur 


GU659691 


Dryas iulia modesta 


MAL-02696 


9-Jul-02 


1 8.608 N 89.846 W 


Mexico: Campeche, EjIdo Conhuas: camino a la zona 
arqueologica Nadzcaan 


GU659594 


Dryas iulia modesta 


MAL-02664 


23-Feb-99 


18.253 N 89.451 W 


Mexico: Campeche, Ejido Narciso Mendoza 


GU659620 


Dryas iulia modesta 


MAL-02664 


23-Feb-99 


18.253 N 89.451 W 


Mexico: Campeche, EjIdo Narciso Mendoza 


GU659622 


Dryas iulia modesta 


MAL-02663 


3-Jan-93 


21.204 N 86.713 W 


Mexico: Quintana Roo, Isia Mujeres 


GU659619 


Dryas iulia modesta 


MAL-02659 


3-NOV-04 


19.157 N 87.544 W 


Mexico: Quintana Roo, Reserva de la Biosfera Sian ka'an: 
Camino a Tampalam 


GU659623 


Dryas iulia modesta 


MAL-02698 


ll-Mar-04 


19.722 N 87.812 W 


Mexico: Quintana Roo, Reserva de la Biosfera Sian ka'an: 
Estacion Santa Teresa 


GU659588 


Dryas iulia modesta 


MAL-02662 


14-Aug-95 


18.367 N 88.585 W 


Mexico: Quintana Roo, Sabldos 


GU659626 


Dryas iulia modesta 


MAL-02665 


20-Jan-07 


21.344 N 87.63 W 


Mexico: Yucatan, TIzlmIn, La Florida 


GU659621 


Dryas iulia modesta 


YB-BCn550 






Panama 


HM416486 


Dryas iulia modesta 


MLA. DL-uo-1033 


1 1-Feb-08 


/./ol 1 N yo./ /U3 b 


mailand: Krabi, Kon rni rni Don 


I A nco c ^ 

KJ4963D2 


Dryas iulia modesta 


MCZ: PBG-08-D001 


l-Jun-08 


7.8834 N 98.3956 E 


Thailand: Phuket, Phuket City (Ampuh Mueng), 
Phuket Butterfly Garden 


KJ496359 


Dryas iulia modesta 


MCZ: PBG-08-D002 


l-Jun-08 


7.8834 N 98.3956 E 


Thailand: Phuket, Phuket City (Ampuh Mueng), 
rnuKGT DUiTerTiy uaruen 


KJ496360 


Dryas iulia modesta 


MCZ: PBG-08-D004 


l-Jun-08 


7.8834 N 98.3956 E 


Thailand: Phuket, Phuket City (Ampuh Mueng), 
Phuket Butterfly Garden 


KJ496361 


Dryas iulia modesta 


MCZ: PBG-08-D005 


l-Jun-08 


7.8834 N 98.3956 E 


Thailand: Phuket, Phuket City (Ampuh Mueng), 
Phuket Butterfly Garden 


KJ496362 


Dryas iulia modesta 


MCZ: PBG-08-D006 


l-Jun-08 


7.8834 N 98.3956 E 


Thailand: Phuket, Phuket City (Ampuh Mueng), 
Phuket Butterfly Garden 


KJ496363 


Dryas iulia modesta 


MCZ: PBG-08-D007 


l-Jun-08 


7.8834 N 98.3956 E 


Thailand: Phuket, Phuket City (Ampuh Mueng), 
Phuket Butterfly Garden 


KJ496364 
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Table 1. Cont. 





Species 


LocBtioni 
Voucher # 


Coll. Date 


Coordinates 


Collection Locality 


GenBank 


Dryas iuHo modesta 


MCZ: PBG-08-D008 


l-Jun-08 


7.SS34 N 98.3956 E 


Thailand: Phuket, Phuket City (Ampuh Mueng), 
Phuket Butterfly Garden 


KJ496365 


Dryas iulia modesta 


MCZ: PBG-08-D009 


l-Jun-08 


7.SS34 N 98.3956 E 


Thailand: Phuket, Phuket City (Ampuh Mueng), 
Phuket Butterfly Garden 


KJ496366 


Dryas iulia modesta 


MCZ: PBG-08-D010 


l-Jun-08 


7.SS34 N 98.3956 E 


Thailand: Phuket, Phuket City (Ampuh Mueng), 
Phuket Butterfly Garden 


KJ496367 


Dryas iulia modesta 


MCZ: LD-08-A15 


1 -Apr-OS 


9.5 N 1 00 E 


Thailand: Surat Thani, Koh Samui 


KJ496353 


Dryas iulia modesta 


MCZ: LD-08-A17 


l-Apr-08 


9.5 N 1 00 E 


Thailand: Surat Thani, Koh Samui 


KJ496354 


Dryas iulia modesta 


MCZ: LD-08-A19 


1 -Apr-OS 


9.5 N 1 00 E 


Thailand: Surat Thani, Koh Samui 


KJ496355 


Dryas iulia modesta 


MCZ: LD-08-A20 


1 -Apr-OS 


9.5 N 1 00 E 


Thailand: Surat Thani, Koh Samui 


KJ496356 


Dryas iulia modesta 


MCZ: LD-08-A21 


1 -Apr-OS 


9.5 N 1 00 E 


Thailand: Surat Thani, Koh Samui 


KJ496357 


Dryas iulia modesta 


MCZ: LD-08-A22 


1 -Apr-OS 


9.5 N 1 00 E 


Thailand: Surat Thani, Koh Samui 


KJ496358 


Dryadula phaetusa 


DHJ: 05-SRNP-31045 






Costa Rica 


GU157119 



GenBank accession numbe rs beginning with KJ correspond to novel sequences generated in this study; all other accession numbers represent sequences downloaded 
from GenBank for inclusion in the analysis. Subspecies identifications of Central American specimens are inferred based on their collection locality. Voucher locations: 
MCZ = Museum of Comparative Zoology, Harvard University; DHJ = D. H. Janzen collection. Some voucher locations were not specified on GenBank. 
doi:l 0.1 371 /journal.pone.Ol 04076.t001 



GAAACAGCTATGACTAAACTTCAGGGTGACCAAAAAAT- 

CA-3'). These primer sequences include the original primers of 
Folmer et al. [21] to which M13 tails (indicated in bold) had been 
concatenated on the 5' end [22]. Addition of these tails to the 
primers increases PGR success, particularly on specimens with 
degraded DNA [23]. PGR products were visualized on agarose 
gels before being sent to Genewiz (genewiz.com) for PGR clean-up 
and bidirectional sequencing. The primer "tails" M13F and 
M13R were used as sequencing primers [22]. We sequenced the 
GOI barcode from 18 D. iulia specimens, constituting all Thai 
and Dominican samples in our dataset. We added 28 additional D. 
iulia barcode sequences from Gosta Rica, Mexico, and Panama 
and a sequence from the outgroup Dryadula phaetusa (Nympha- 
lidae: Heliconiinae) to the genetic dataset. We included all Dijas 
iulia sequences in GenBank that completely overlapped with the 
barcoding fragment that we sequenced; longer sequences were 
trimmed so that each sample included exactly 658 bp. The two 
Dominican specimens represent the subspecies D. iulia hispaniola. 
AH other sequences, including those from Thailand, are of Z). iulia 
modesta. Sequences from the other eleven D. iulia subspecies — 
which are mosdy Caribbean island endemics — were unavailable. 
Sequences were viewed, assembled, aligned, and trimmed with 
Geneious [24]; alignments were performed within Geneious using 
MUSCLE [25]. The sequence alignment is provided as a nexus 
file in Appendix SI. Protocols were adopted from dnabarcodin- 
gl01.org, developed by the DNA Learning Center, Cold Spring 
Harbor Laboratory. 

Phylogenetic and distance analyses 

The most parsimonious haplotype network of D. iulia was 
determined with TCS 1.2 with a 95% connection limit [26], and 
redrawn using the Pie Graph Tool in Adobe Illustrator CS6 
(adobe.com). The program jModelTest 2.1.4 [27] was used to 
select the GTR-I-I-I-G model of sequence evolution using the AIC 
criterion, but we implemented the GTR-I-G model to avoid 
overparameterizing the data. A maximum likelihood analysis and 
an ML rapid bootstrap analysis were performed with RAxML 
7.6.3 [28] on the CIPRES Science Gateway (phylo.org) [29]. 
Bootstrapping was stopped automatically using the majority rule 



criterion under a GTR-I-G model. Bayesian phylogenetic analyses 
were performed with MrBayes 3.2.2 [30] on the GIPRES Science 
Gateway. Four Markov chains, one cold and three heated, were 
run simultaneously for 10 million generations. Trees were sampled 
every 1000 generations, and the first 25% of sampled trees were 
discarded as burn-in before calculating a consensus tree. Changes 
in the posterior probabilities of 20 nodes were plotted over the 
generations of the analysis with the program Ai'e We There Yet? 
[3 1] in order to confirm that the chains had probably converged. 
To assess parsimony support for relationships among taxa, 1000 
bootstrap replicates were run in TNT 1 . 1 using standard 
bootstrapping with replacement after "Max. trees" was reset to 
10,000 [32]. TaxonDNA 1.0 (taxondna.sourceforge.net) [33] was 
used to calculate uncorrected p-distances between barcode 
sequences. 

Results 

DNA barcode sequences were identical among all specimens 
from Thailand: the single wild-caught Dryas iulia on Koh Phi Phi 
Don, the six wild-caught specimens on Koh Samui, and all nine 
specimens obtained from the Phuket Butterfly Garden (Fig. 2a), 
This 658 bp haplotype was not shared with any specimens caught 
in the New World, but was most similar to a Costa Rican 
specimen (1 bp difference). The tree topologies obtained from 
Bayesian, maximum likelihood, and parsimony methods were 
similar and had universally poor branch support, as one might 
expect of a phylogeny based on a single gene sampled within a 
single species (Fig. 2b). Each of the two D. iulia hispaniola 
specimens sampled from the Dominican Republic had a unique 
haplotype; both were notably distinct from the other sampled 
haplotypes. Several haplotypes were found in both Mexico and 
Costa Rica, demonstrating genetic diversity within the subspecies 
D. iulia modesta. 

Discussion 

It is likely that wild populations of Dryas iulia in Thailand 
originated from livestock at a butterfly farm, possibly from 
individuals that were intentionally released alive. One such farm. 



PLOS ONE I www.plosone.org 



4 



August 2014 | Volume 9 | Issue 8 | e104076 



DNA Barcodes Identify the Origin of an Alien Species 



the Phuket Butterfly Garden (PBG), may be the source of the wild 
population, but without exhaustive sampling, we cannot rule out 
other such farms as potential sources of the naturalized wild Thai 
population. There are at least two alternative scenarios consistent 
with our results. It is possible that the D. iulia livestock at PBG was 
obtained from the same source as a second, unknown source that 
was responsible for the introduction — perhaps another butterfly 
farm in Southeast Asia unknown to us. Alternatively, D. iulia 
could have been introduced into the wild in Thailand where they 
became established and subsequently collected by PBG for 
propagation and sale. We consider both of these alternatives 
unlikely. We know of no other butterfly houses in Thailand that 
stock non-Asian species, including Nong Nooch Tropical Garden, 
Pattaya (nongnoochgarden.com), Siam Insect Zoo, Chiang Mai 
(malaeng.com), Bangkok Butterfly Garden and Insectarium, and 
Bai Orchid and Butterfly Garden in Chiang Mai. Access to import 
documentation or knowledgeable PBG staff members could 
confirm or refute the second possibility. 

Identical sequences between wild-caught specimens and those 
from PBG are not due to lack of genetic diversity in the species or 
subspecies. We included all publicly available, homologous D. 
iulia barcode sequences in our dataset (which happened to all be 
from the subspecies D. iulia modesta), and the 30 sequences from 
non-Thai samples constitute 1 7 distinct haplotypes, demonstrating 
some degree of genetic variability within the species as a whole. 
The lack of genetic diversity within Thai D. iulia is consistent with 
a genetic bottieneck caused by a small founding population. This 
might have happened if a small number of individuals was 
imported to Thailand and used to found a colony at a butterfly 
house that eventually became inbred [34] . 

It is unclear how this novel introduction wiU affect wild 
populations of other organisms. The species has been observed 
feeding on Passiflora foetida, which is an invasive plant in 
Thailand, and the butterfly might therefore be a boon for 
biological control of this weed. However, herbivory by D. iulia 
might suppress populations of other species through consumptive 
competition. This vine also provides fodder for the native butterfly 
species Cethosia cyane and Vindula erota, as well as the alien 
species Acraea terpsicore [ = A. violae) [13], which could be 
adversely afferted. 

Naturalization of this exotic species in Thailand may not have 
dire ecological consequences. The Monarch butterfly, Danaus 
plexippus, dispersed across the Pacific Ocean from the Americas to 
Australia in the IQ*^ century [35,36]. This relatively recent 
addition to Australia's biota does not seem to suppress populations 
of native species, particularly since the larvae feed on introduced 
plant species including Asclepias curassavica and Gomphocarpus 
fruticosus (-Asclepias fruticosa) [5]. Introduced insect species 
occasionally increase their host breadth to include plant species 
native to the area of introduction [37], andZ). iulia might impact 
native vegetation if this occurs. Observations of the species are 
currently confined to peninsular Thailand and Malaysia. Wild D. 
iulia was first recorded in Asia only seven years ago, and the 
species may stiU be expanding its range. Continued live butterfly 
release at weddings and religious ceremonies may be fortifying 
wild populations and aiding range expansion. 

There are several ecological dangers associated with butterfly 
houses. Most of these facilities do not breed butterflies for display. 
Instead, they are sent shipments of live pupae from butterfly farms 
by express maU. Many of these shipments cross international 
borders, as the majority of butterfly farms are in tropical countries 
and many butterfly houses are in temperate areas [18]. A relatively 
small number of butterfly farms supply pupae for most of the 
world's butterfly houses, with large numbers of butterfly farms in 



Costa Rica, Malaysia, and the Philippines. This translocation of 
fivestock opens the possibility that exotic species could escape into 
areas where they are not native, thereby introducing novel and 

potentially pestiferous lepidopteran species into natural ecosys- 
tems, as seems to have happened in Thailand. If an escapee is from 
a species found locally, interbreeding between introduced and 
native genotypes could disrupt locally co-adapted gene complexes. 
This insidious "biopoUution" of a gene pool could be harmful to 
species with separate populations that are locally adapted to 
different conditions. Even if butterflies remain contained within 
the facilities designed to house them, lepidopteran parasites and 
pathogens harbored by the hving, translocated pupae are smaller 
and not easily detected. Escape of these butterfly enemies into the 
wild could have profoundly negative consequences on local 
butterfly populations [18]. However, it is possible that these risks 
can lie offset to some degree by the potential for butterffy houses to 
educate the public about basic biology and the importance of wild 
insects and their habitats, which are threatened around the world 
[18]. 

In early 2014, PBG's website showed photographs of at least 
three different couples in wedding garb releasing butterflies from a 
cage, and at least one, live D. iulia specimen can be seen in each 
photograph (phuketbutterfly.com/wedding.php, accessed March 
2014). In addition to release at weddings, thousands of butterffies 
are released annually into the Khao Phra Thaeo Wildlife 
Conservation Area in Phuket in a release ceremony orchestrated 
in part by PBG [16,17]. The Phuket Gazette, a local newspaper, 
has recorded videos of these events in which release of D. iulia can 
be observed [16,17]. In many parts of Asia, captive animals have 
been released into the wild for over 1,000 years as part of Buddhist 
rituals aimed at cultivating compassion for living beings [38]. In 
recent decades, exotic species are readily available in live animal 
markets in Asia either as pets or food. Release of tlu-sc non-native 
species has led to their establishment as invasive species in some 
areas [38-40]. For example, the American turtie Trachemys 
scripta, which is sold as food and frequentiy released into the wild, 
is now the most common turtie in every river in Taiwan [41]. 
There are several initiatives to educate Buddhist monks and laity 
about the ecological dangers of animal release [38,40,41]. 

Whereas butterfly houses offer the advantages of conservation 
awareness and general education about the importance of 
biodiversity, there are few, if any, positive environmental aspects 
of intentional butterffy release. Species introduction, biopoUution 
of natural gene pools, and introduction of novel butterffy enemies 
are all far more likely when fecund, living butterffies are 
intentionally released into the wild [18,19]. For these and other 
reasons, several authors have called for a ban on the release of 
butterflies at weddings [42-45]. Within the United States, USDA- 
APHIS releases specific guidelines regarding the butterffy species 
that can be legally r(;l(^ased in each stat(; [46] in order to reduce 
the likelihood of negative ramifications of live butterfly release. We 
concur with other authors [42-45] that the release of live 
butterflies at social or cultural events should be banned; the 
short-hved benefits do not justify the threats of long-term damage. 

Conclusions 

Our analyses suggest that the Neotropical butterfly Dryas iulia 
was introduced to Thailand by th(; Phuket Butterfly Garden 
(PBG), which breeds the species for live release at weddings and 
other public events. Most wild D. iulia locality records in Thailand 
and Malaysia are on tropical islands that are frequently the site of 
destination weddings: Phuket, Phi Phi Don, Samui, and Tioman. 
It is likely that PBG obtained hvestock from a butterfly farm in 
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Figure 2. Relationships among Dryas iulia showing that wild-caught. Thai specimens have identical barcodes with specimens from 
the Phuket Butterfly Garden and no other samples collected in Central America and the Caribbean, a) Most parsimonious haplotype 
network of D. iulia constructed with 95% connection limit. The size of each circle is proportional to the number of specimens sharing that haplotype; 
the smallest circles represent a single haplotype and the largest, sixteen. The colors of the pie charts indicate proportional representation of the 
provenance of samples with that haplotype. Separation by a line indicates a single base pair difference between haplotypes; crosses represent 
haplotypes that would be 1 bp different than adjacent haplotypes, but were not sampled in this study, b) maximum likelihood bootstrap consensus 
tree of COI haplotypes from D. iulia and one outgroup. Codes refer to GenBank Accession Numbers and colors denote provenance of specimen 
collection. Numbers near selected nodes indicate refer to the following branch support values (maximum likelihood bootstrap support, Bayes 
posterior probability, parsimony bootstrap support, respectively): 0 = 40, -, -; 0 = 27, 0.62, -; © = 20, -, -; 0 = 27, 0.75, -; 0=31, 0.98, -; 0 = 26, 
0.69, -;©=100, 1, 100; ©=100, 1, 100. 
doi:1 0.1 371 /journal.pone.01 04076.g002 



Costa Rica (as evidenced by similarity of barcode sequences), and 
subsequent inbreeding at tlie PBG expunged genetic variation, if 
there was any in the founding population. Released specimens 
bred in tlie wild and began using Passiflora foetida as a larval host 
plant. The distribution of the species in Thailand currently 
encompasses thousands of square kilometers, and eradication 
efforts are unlikely to be successful, particularly since P. foetida is a 
common, invasive species, making it difficult to find all possible 



larval host plants for control purposes. To strengthen our 
conclusions regarding the provenance of the Thai stock, future 
studies might include more markers and obtain samples from the 
species' entire native range, which includes most islands of the 
Caribbean as well as northern South America. We suggest that 
Thai authorities prohibit the intentional release of live butterflies 
for commercial purposes and social functions, and regulate the 
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importation of live animals for non-scientific purposes, including 
insects, to prevent similar introductions in the future. 

Supporting Information 

Appendix SI Nexus file including the DNA sequence alignment 
analyzed in this paper. Sequence names refer to GenBank 
accession numbers followed by the two-letter abbreviation for 
the country where the specimen was collected. 
(NEX) 
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